Heterologous production and characterization of a pyomelanin of Antarctic Pseudomonas sp. ANT_H4: a metabolite protecting against UV and free radicals, interacting with iron from minerals and exhibiting priming properties toward plant hairy roots

Background Antarctica has one of the most extreme environments in the world. This region is inhabited by specifically adapted microorganisms that produce various unique secondary metabolites (e.g. pigments) enabling their survival under the harsh environmental conditions. It was already shown that these natural, biologically active molecules may find application in various fields of biotechnology. Results In this study, a cold-active brown-pigment-producing Pseudomonas sp. ANT_H4 strain was characterized. In-depth genomic analysis combined with the application of a fosmid expression system revealed two different pathways of melanin-like compounds biosynthesis by the ANT_H4 strain. The chromatographic behavior and Fourier-transform infrared spectroscopic analyses allowed for the identification of the extracted melanin-like compound as a pyomelanin. Furthermore, optimization of the production and thorough functional analyses of the pyomelanin were performed to test its usability in biotechnology. It was confirmed that ANT_H4-derived pyomelanin increases the sun protection factor, enables scavenging of free radicals, and interacts with the iron from minerals. Moreover, it was shown for the first time that pyomelanin exhibits priming properties toward Calendula officinalis hairy roots in in vitro cultures. Conclusions Results of the study indicate the significant biotechnological potential of ANT_H4-derived pyomelanin and open opportunities for future applications. Taking into account protective features of analyzed pyomelanin it may be potentially used in medical biotechnology and cosmetology. Especially interesting was showing that pyomelanin exhibits priming properties toward hairy roots, which creates a perspective for its usage for the development of novel and sustainable agrotechnical solutions. Supplementary Information The online version contains supplementary material available at 10.1186/s12934-022-01990-3.


Background
Antarctica is one of the most extreme regions on Earth, characterized by very low temperatures, strong ultraviolet (UV) radiation, and sunlight and nutrient deficiency [1]. However, even harsh polar environments are inhabited by highly specialized microorganisms that have developed a number of adaptive traits necessary for their survival [2]. These include the synthesis of various secondary metabolites, e.g., biosurfactants and siderophores, which help to increase the availability of nutrients and minerals [3] or the synthesis of pigments (carotenoids and melanins) responsible for protection against UV radiation [4,5].
An important group of microorganisms often isolated from the Antarctic regions is bacteria of the genus Pseudomonas [6]. These rod-shaped, Gram-negative bacteria are valued in biotechnology due to their high metabolic potential for the production of several valuable secondary metabolites [7]. Among many known compounds, pseudomonads produce bioactive pigments classified as phenazines (pyocyanin and pyorubrin), siderophores (pyoverdine, fluorescein), and melanins (pyomelanin) [8].
While much is known about the basic functions of the above-mentioned secondary metabolites, the potential of melanins is still not fully exploited.
Melanins are brown-black polymers produced by organisms mainly to protect against UV radiation [9]. These bioactive polymers are classified according to the chemical precursors used in their biosynthesis, being eumelanin, pheomelanin, neuromelanin, allomelanin, and pyomelanin, the latter of which is also classified as a form of allomelanin [10,11]. Pyomelanin originates from the catabolism of L-tyrosine or L-phenylalanine which, due to further transformation by 4-hydroxyphenylpyruvate dioxygenase (Hpd), forms homogentisic acid (HGA). HGA is the substrate of homogentisate 1,2-dioxygenase (HmgA). However, overexpression of Hpd or repression of HmgA leads to accumulation of HGA and causes its secretion from the cell where it undergoes oxidation to benzoquinoneacetic acid and self-polymerization into pyomelanin [10,12]. Moreover, due to the heterogenous nature of the polymer [13] and under certain circumstances, additional substrates may appear during its formation as a consequence of, e.g., quinone oxidoreductases activity [14][15][16].
The chemical structure of pyomelanin is responsible for numerous biological properties and functions such as protection from light and oxidative stress, energy transduction as well as chelation and reduction of metal ions, e.g., Fe 3+ [10,17]. It has been suggested that pyomelanins may contribute to increasing the availability of certain elements in minerals for the microorganisms [9,10,17]. They can, therefore, play a supporting role for bacterial siderophores through increasing iron bioavailability [18]. On the other hand, since melanins produced by certain plant pathogens are virulence factors [19], these compounds, as elicitors, can stimulate defense responses in plants [20]. Checking whether a certain compound can be used as a potential elicitation agent very often involves examining its effect on the production of general and specialized metabolites in native plants or their in vitro cultures. The results are often manifested in the modified proportions of both types of metabolites [21]. In summary, these data confirm that melanins have considerable biotechnological potential.
In this study, the identification and functional characterization of pyomelanin originating from the Antarctic, psychrotolerant Pseudomonas sp. ANT_H4 was performed. The subsequent optimization of the pyomelanin production allowed revealing biological functions of this pigment as well as its potential for being used in environmental biotechnology.

General physiological characterization of Pseudomonas sp. ANT_H4
Pseudomonas sp. ANT_H4 shows growth at temperatures from 4 to 37 °C (optimum 18 °C), tolerates a pH of 6-10, and can survive in NaCl concentrations of up to 3% [6]. It is also worth mentioning that ANT_H4 exhibits many features considered to be environmental adaptation strategies: ANT_H4 is able to produce siderophores in amounts of ~ 0.7 mmol/l as well as biosurfactants, which are able to lower the interfacial tension of a medium up to ~ 55% (from 55 (± 2) to 30 (2 ±) mN/m). An important feature distinguishing ANT_H4 from most of the other Antarctic Pseudomonas strains studied is its ability to sequester a brown-black pigment into the medium, in amounts of ~ 0.01 mg/mL.

Overproduction of melanin-like compounds from fosmid clones
From among more than 3000 E. coli MPO554 fosmid clones carrying ~ 45 kb genome fragments of Pseudomonas sp. ANT_H4 genome, we selected those that produced brown-black pigment into the medium. The expression medium supplemented with L-tyrosine allowed fosmid clones to efficient production of a brown-black pigment in amount of 0.25 mg/ mL (Fig. 1). The obtained metabolite was compared chromatographically to that obtained from the native strain ANT_H4 and revealed identical nature. This is the first time that an overproduction of a melanin-like compound has been achieved through a fosmid expression system. The single-gene heterologous expression systems constructed by other scientists enabled to receive ~ 0.213 mg/mL from E. coli JM109 [22] and ~ 0.315 mg/mL from E. coil KSYH [23] using the 4-hydroxyphenylpiruvate dioxygenase gene from P. aeruginosa and Ralstonia pickettii, respectively.

Identification of a melanin-like compound produced by Pseudomonas sp. ANT_H4
As mentioned above Pseudomonas sp. ANT_H4 showed the unique ability to sequester a brown-black pigment into the medium. Based on the literature review, we concluded that this may be a melanin-like compound [24][25][26]. For its identification, HPLC-DAD and FT-IR analyses were performed.
Analysis by HPLC using a DAD detector allowed for the determination of three basic characteristics of the compound. First, the comparatively high polarity due to poly-quinone pyomelanin structure [10], as evidenced by a relatively short retention time (R t = 4.2 min, by Fig. 2A) in the water/acetonitrile gradient. Second, the maximum UV absorbance of the compound at a wavelength of λ = 230 nm (Fig. 2B), which is a characteristic feature of compounds from the melanin group [27]. Three, a very broad peak typical for polymers [28].
In addition, the FT-IR spectroscopy of a sample of the ANT_H4-derived melanin-like metabolite was performed. The analysis of the IR spectra revealed several melanin-characteristic signals [29][30][31], i.e., a broad band at 3434 cm −1 assigned to the presence of an OH group (tensile vibration), a band at 1728 cm −1 corresponding to vibrational excitation of carbonyl groups C = O or of COOH groups, a band at 1634 cm −1 may come from the ionized COO − groups, at 1403 cm −1 indicating vibration from ionized COO − groups, and at 1241 cm −1 that may be a signal from the combination of C-O bonds (Fig. 3).
Altogether, the performed HPLC and FTIR analyses confirmed the presence of a melanin-like compound classified as heterogenous pyomelanin [32], which is produced by some Pseudomonas species [24][25][26]. Overproduction of melanin-like compound in E. coli MPO554 carrying fosmid with appropriate ANT_H4 genome fragment. The left flask shows the uninoculated LB medium, while the right flask presents the ANT_H4 culture cultivated for 120 h using the same medium Fig. 2 Analytical characterization of the ANT_H4-derived melanin-like compound. a reversed-phase C18 HPLC chromatogram, revealing a broad peak with a retention time Rt = 4.2 min, and b UV/Vis spectrum with absorption maxima at λ = 230 nm Genomic characterization of Pseudomonas sp. ANT_H4tracking of genes responsible for pyomelanin production Sequencing of the Pseudomonas sp. ANT_H4 genome using the Illumina MiSeq platform generated 1,030,406 paired reads and 319 Mbp of sequence information. As a result of the assembly of the Pseudomonas sp. ANT_H4 genome, 89 contigs of a total length of 6,131,973 bp were obtained. Initially, the genome sequence was automatically annotated using RAST on the PATRIC 3.6.2 web service and its general features for the ANT_H4 strain are presented in Table 1.
To identify genes directly responsible for the production of pyomelanin, functional analysis using a fosmidbased expression system was applied [41]. Four fosmid clones displaying the capacity of synthesizing pyomelanin were further sequenced. The analysis of these sequencing results revealed a common core set of 27 genes (Fig. 4, Table 2). Amongst the identified enzymes, we found one gene function that might be directly involved in the melanin-overproducing phenotype: quinone oxidoreductase (EC 1.6.5.5). The gene encoding for this enzyme was amplified by PCR, subcloned into plasmid pBluescript II SK, and transformed into E. coli TG1 and MPO554. However, the pyomelanin production phenotype was not observed in any of the tested hosts, which may indicate that either the introduced gene is inactive, not involved in pyomelanin production, or not sufficient for pyomelanin overproduction. Taking into account that all fosmids in pyomelanin-producing strains contained a set of common genes we assume that this is not a single-gene-generated phenotype, but rather a consequence of the action of several genetic determinants creating a unique cluster.
We further explored the presence and co-location of genes encoding for proteins homologous to proteins encoded by the core of fosmids giving the desired phenotype. After a search of the NCBI nr database, we observed that there were 64 contigs carrying at least 10 out of 23 homologous genes and 34 of them carried these genes within a 50-kb-range fragment of DNA (Fig. 4). Besides the ANT_H4 contig (accession number VOBG01000010), we identified the same gene cluster in another strain isolated from the same location -Pseudomonas sp. ANT_H12B (VOBH01000017). Other strains encoded at most 15 homologous proteins. It is also worth noting that within the set of genomic fragments, only P. lurida PGSB 7828 (D010000004) encoded for a quinone oxidoreductase homolog.

Evaluating the biotechnological potential of the identified pyomelanin Sun-protection and free radical scavenging activity
One of the basic functions of the compounds from the melanin group is protection against UV radiation. The ability of ANT_H4-derived pyomelanin to absorb UV light contributes to direct protection from sunlight. As a non-cytotoxic compound [15,25,42], pyomelanin could be an excellent additive to enhance the power of sunscreens. In order to test that hypothesis, the sun protection factor (SPF) of three commercial sunscreens was determined empirically: The stated SPFs were 5, 15, and 25, while the results of the analyses revealed values of 5.09 ± 0.03, 13.41 ± 0.35, and 25.10 ± 0.13, respectively. Addition of pyomelanin to each of these sunscreens resulted in an increase of SPF to: (i) 7.79 ± 0.04, 15.88 ± 0.43, 27.48 ± 0.12 when pyomelanin in final concentration of 0.1 mg/mL was added, and (ii) 6.41 ± 0.14, 14.65 ± 0.19, 25.35 ± 0.17 when pyomelanin in final concentration of 0.01 mg/mL was added [25,43,44]. The effect of ANT_H4-derived pyomelanin is consistent with the results of other researchers testing the SPF changes of sunscreen [25,42,45]. However, currently, only two melanin-like compounds are used in cosmetics; catechol melanin from black olive fruits (Laboratoires Biocyte, France), and "bio-melanin" in anti-UV lotions (Baar Products, PA, USA). So far, there is no mention of the use of pyomelanin, mainly due to the lack of optimization of its production efficiency [42].
UV radiation (especially UVA) reaching the skin leads to the generation of free radicals that oxidize proteins, lipids, and nucleic acids [46]. In the case of pyomelanin, the ability to scavenge free radicals is its primary feature [24,47]. In order to test the antioxidant capacity of the ANT_H4-derived pyomelanin, the free-radical scavenging effect was tested on a 0.1 mM DPPH solution (used as a free radical). The IC 50 value (antioxidant concentration required for quenching 50% of the initial DPPH) was established at a concentration of ~ 0.75 mg/mL pyomelanin, which quenched maximally 60% of the DPPH at a concentration of 2 mg/mL (Fig. 5). This is in line with the results of previous studies [15,17,45].
It is worth mentioning that the antioxidant activity significantly increases the applicability of the tested compound in cosmetology. This feature is highly desirable in, e.g., creams or lotions that reduce the negative impact of free radicals on skin aging.

Iron binding properties
One of the features of pyomelanin is its ability to interact with iron-containing particles. The presence of quinones in the polymer structure, as well as the specific structure resembling siderophores in certain Fig. 4 A comparison of the fosmid-encoded core proteome with other bacterial genomes. The latter were selected based on the presence of at least 10 homologous proteins encoded within a 50-kb-range of DNA sequence. From these, only homologous regions were extracted based on coordinates present in brackets below the accession number of a particular contig. They were aligned with a clinker tool, where each protein coding sequence is represented as a rectangle arrow. Homologous proteins have the same color and proteins from adjacent contigs are connected with a block of a lighter equivalent of their color. The genus of the host bacterial strain is color-coded next to accession numbers. Fragments originating from the genome of ANT_H4 were marked with a dashed, red frame. Additionally, for the sake of figure clarity, only selected proteins from the fosmid-encoded core proteome were labeled with an accession number or locus_tag configurations, enables both the binding and reduction (Fe 3+ to Fe 2+ ) of iron [10,48]. For example, pyomelanin produced by Legionella pneumophila conferred a ferric reductase activity and thus played an important role in iron uptake [18,30,49]. This is a feature that leads to an increase in iron bioavailability and thus positively impacts bacterial fitness. The iron-binding properties of ANT_H4-derived pyomelanin were performed against the CAS-HDTMA-iron complex, in which iron is relatively weakly bonded [50]. The iron chelation capacity of pyomelanin was measured spectrophotometrically at OD 630 . The results were similar to deferoxamine mesylate siderophore at a concentration of 0.025 mM and reached its maximum at a concentration of 2 mg/mL. Therefore, it suggests the participation of the tested compound in increasing the availability of iron in the environment, especially as the structure of the pyomelanin and its iron binding mechanism are similar to humic acids, commonly used in bio-fertilizers for plants [51][52][53][54].
Many iron chelators also exhibit antioxidant properties [55]. It was shown that the iron-binding and free-radical scavenging properties of ANT_H4-derived pyomelanin correlated and reached their maximum at a concentration of 2 mg/mL of pyomelanin (Fig. 5). This result is very important according to the iron-reduction properties of the investigated compound, which can lead to Fenton's reaction, responsible for the generation of free radicals [54,56].

Bioconsolidation with ANT_H4-derived pyomelanin
Considering the pyomelanin iron-binding properties, its influence on minerals was tested. Bioconsolidation studies were conducted with the use of pyrite and hematite. Their basic composition is FeS 2 and Fe 2 O 3 , i.e., the iron is in a reduced (Fe 2+ ) and oxidized (Fe 3+ ) form, respectively.
The batching experiment revealed increased iron retention in the presence of ANT_H4-derived pyomelanin. Both tested concentrations, i.e., 0.1 and 0.01 mg/mL significantly reduced the loss of iron after 1, 24, and 48 h of batching minerals. The results after 48 h were similar for both pyomelanin concentrations and showed over 80-fold increased hematite iron retention compared to the control. In the case of pyrite, an increased (over 60-fold) retention was observed with the use of a higher (0.1 mg/mL) concentration of pyomelanin (Table 3). The results are very promising in light of ANT_H4-derived pyomelanin applicability. Due to the specific chemical structure of pyomelanin, it participates in iron redox processes, but we also hypothesize that the polymer may adhere to the surface of minerals and stabilize them by creating a protective layer.

Hairy roots elicitation
Melanin-like compounds are often identified among many fungal plant pathogens as their virulence factors protecting non-specifically against host defense mechanisms. The pathogen-free molecules themselves are not harmful to plants, however, they are recognized by them as a potential threat [47]. This initiates a cascade of biochemical processes leading to modifications of biosynthetic pathways resulting in different proportions of plant general and specialized metabolites, and thus to a preventive enhancement of the systemic defense response [21]. Sterols are considered to be plant general metabolites as they build cell membranes and are in charge of maintaining the fluidity and permeability of these structures. Phytosterols also play a role in adaptation to many stress (mostly abiotic) conditions [57]. Pentacyclic triterpenoids are plant specialized metabolites that take part in defense strategies against herbivores and pathogens [58]. During stressful conditions as well as upon elicitation, biosynthesis pathways of sterols and pentacyclic triterpenoids may be competitive due to the common precursor of both groups of compounds, i.e., squalene [21]. However, some factors may stimulate both pathways in parallel, showing evidence that plants' defense strategies are combined and factor-dependent.
The sterol fraction obtained from diethyl ether extracts was composed of a group of typical sterols, namely campesterol, cholesterol, isofucosterol, sitosterol, and the predominating stigmasterol. One steroid ketone, tremulone, and one biosynthetic precursor, 24-methylenecycloartanol, were also identified in this fraction (Additional file 1: Table S1). The presence of pyomelanin in the culture medium caused a significant increase in steroid biosynthesis; the total content of these compounds in pyomelanin-elicited hairy roots was approximately 58% and 24% higher than in the control roots after 7 and 14 days, respectively (Fig. 6, Additional file 1: Table S2).
Four sterols, i.e., cholesterol, campesterol, sitosterol, and stigmasterol, were also found in glycosides of methanol extracts subjected to hydrolysis (Fig. 7, Additional file 1: Table S3). The contents of sterol glycosides were similar in the control and elicited samples. After seven days of incubation, the number of sterols slightly increased (13%) in the elicited samples, however, after 14 days it decreased (22%) in comparison with control. The modifications in the ratio of conjugated to free forms of sterols are often considered a major key to modulating the membrane order and membrane biophysical properties. The changes in the content of sterol glycosides are often observed during plant stress reactions due to their influence on such plant membrane properties as fluidity and permeability. These changes can also influence the activity of some membrane-bound enzymes as well as proton pumps, which can lead to further modifications of various metabolic pathways [59].
A notable effect of pyomelanin elicitation was observed in the content of free oleanolic acid in hairy roots tissueit increased more than 18-fold after 7 days and more than 90-fold after 14 days in the elicited samples (Fig. 8, Additional file 1: Table S4). Simultaneously, the content of the oleanolic acid precursor, triterpenoid alcohol β-amyrin, increased almost tenfold after 14 days of elicitation (Additional file 1: Table S2). The number of oleanolic acid saponins released into the medium was also higher in pyomelanin-treated samples, but the differences were not statistically significant (Fig. 9A, Additional file 1: Table S5). At the same time, the content of oleanolic acid saponins accumulated in the hairy roots tissue of the elicited samples was around 2.7-and 2.9-fold lower as compared with the control after 7 and 14 days, respectively (Fig. 9B, Additional file 1: Table S6). This might be explained by the increased release to the medium, however, the observed enhanced saponin secretion does not fully compensate for the decrease in their accumulation in root tissue occurring after elicitation. Therefore, it  seems that it is the free oleanolic acid, and not its saponins, that plays the main role in Calendula officinalis Linnaeus hairy roots defense strategy against potential virulence factors. The diverse effects of several elicitors including chitosan, heavy metals, and phytohormones, on the biosynthesis of sterols and triterpenoids in Calendula officinalis L. hairy roots, have been reported previously [59,60]. The most spectacular effect observed so far was the 113fold increase of oleanolic acid saponin release to the medium after elicitation with jasmonic acid [61]. However, such a high boost to the synthesis of free oleanolic acid as demonstrated in the present study (a more than 90-fold increase after elicitation with ANT_H4-derived pyomelanin) is observed for the first time and has never been reported previously. This result indicates that ANT_H4-derived pyomelanin elicitation is particularly effective in the stimulation of biosynthesis of triterpenoids occurring in a free, and not glycosidic, form. Such specificity in the modulation of the triterpenoid biosynthetic pathway seems to be a rather unique feature of the applied elicitor, and it might be of particular importance for the potential enhancement of biosynthesis of other bioactive triterpenoids.
The obtained results point to the possibility of the application of received pyomelanin in agrotechnics, e.g., as a potential plant priming agent. However, the use of this type of elicitor still requires additional studies of native plants as well as on a larger scale.

Conclusions
Nowadays the search for natural, biologically active molecules is highly needed. Bioprospection of microorganisms adapted to extreme environmental conditions offers great opportunities for finding new biotechnological solutions. The selection of the Antarctic Pseudomonas sp. ANT_H4 strain that produce pyomelanin into the environment, as well as establishing its efficient overproduction using the fosmid expression system, made it possible considering the applicability of extracted biopolymer. Besides the basic, confirmed properties of ANT_H4-derived pyomelanin, such as sun and free radical protection and interaction with iron ions, we revealed its priming abilities to in vitro hairy root cultures. These interesting results requires further research on native plants in their natural environment and thus, subsequent studies will be directed into development of novel and sustainable agrotechnical solutions based on pyomelanin.

Bacterial strains, vectors, and culture conditions
Pseudomonas sp. ANT_H4 was obtained from a collection of bacterial cultures that were previously isolated from soil samples taken in 2012 from King George Island (Antarctica; GPS coordinates: 62 09.6010 S, 58 28.4640 W) [6]. E. coli MPO554 and pMPO579 vectors were used for fosmid library construction [41]. For subcloning, E. coli TG1 [62] along with plasmid pBluescript II SK [63] were used.
All strains were grown in LB medium in aerial conditions (3 RCF, Relative Centrifugal Force). Pseudomonas sp. ANT_H4 was grown at 18 °C while E. coli MPO554 and E. coli TG1 were grown at 37 °C. For the fosmid expression, the medium was supplemented with 12.5 µg/ mL chloramphenicol, 800 µg/mL salicylate, 100 µg/mL arabinose and 1 mg/mL L-tyrosine. To obtain an expression of genes cloned in fosmid, it was supplemented with 50 μg/mL X-gal, 120 µg/mL IPTG, and 100 µg/mL ampicillin. When necessary, media were solidified by the addition of 1.5% (w/v) agar.

Hairy roots cultures
Calendula officinalis L. hairy root line CC16 was obtained from cotyledon explants infected with Agrobacterium rhizogenes strain ATCC 15,834 according to a previously described procedure [64]. The transformation was confirmed by GUS (β-glucuronidase) reporter gene system in the histochemical assay and by PCR analysis of the rolC gene. The roots were cultivated in a liquid-modified Murashige and Skoog (MS) medium [65] with halved macronutrient concentrations (½ MS, at 25 °C, in the dark on a rotatory shaker at 1.5 RCF). Subcultures were performed every 3 weeks by transferring 2 cm pieces of the young, branched root to 100 mL of a fresh medium.

Analysis of the production of siderophores
To determine the number of siderophores produced, bacteria were cultivated for 7 days in a GASN medium at 18 °C with rotary shaking set to 3 RCF. The initial optical density at 600 nm (OD 600 ) was 0.1. After 7 days of cultivation, bacteria were centrifuged (4000 RCF for 5 min) and supernatants were added in a 1:1 ratio to the CAS reagent (CAS-HDTMA-chrome azurol S-hexadecyltrimethylammonium bromide) [66]. GASN medium was used as a negative control, while deferoxamine mesylate salt (Sigma-Aldrich, Saint Louis, MI, USA), at a concentration of 0.025 mM, was used as a positive control. All experiments were performed in triplicate. After an hour of incubation, the absorbance at 630 nm was measured using an automated microplate reader.

Analysis of the production of biosurfactants
To determine the presence of biosurfactants, bacterial cultures were cultivated for 7 days in LB medium at 18 °C with rotary shaking set to 3 RCF in two variants: (i) with 1% (w/v) sunflower oil and (ii) without this supplementation. The initial optical density at 600 nm (OD 600 ) was 0.1 and the controls were uninoculated variants. The experiments were performed in triplicate. After 7 days of cultivation, bacteria were centrifuged (4000 RCF for 5 min) and supernatants were tested through the ring method (du Nouy method) using a Kruss Tensiometer K20 (Kruss GmbH, Hamburg, Germany) [67].

Draft genome sequencing
Genomic DNA of the Pseudomonas sp. ANT_H4 was isolated using the cetyltrimethylammonium bromide (CTAB)/lysozyme method [68]. Illumina TruSeq libraries for each strain were constructed following the manufacturer's instructions. The genomic libraries were sequenced on an Illumina MiSeq instrument (using the v3 chemistry kit) (Illumina, San Diego, CA, United States) in the DNA Sequencing and Oligonucleotide Synthesis Laboratory (oligo.pl) at the Institute of Biochemistry and Biophysics, Polish Academy of Sciences, Warsaw. Raw reads were processed using fastp [69] version 0.19.5 with the following flags: -cut_by_quality3 -cut_window_ size 10 -cut_mean_quality 25 -trim_poly_x -poly_x_ min_len 15 -length_required 100. Filtered reads were used for assembly with SPAdes version 3.11.1 with -careful flag.

Bioinformatic analysis
The analyzed bacterial genomes were automatically annotated using RAST [70] on the PATRIC 3.6.8 [71] web service and manually curated. Similarity searches were performed using BLAST programs [72]. The metabolic features were identified with the SEED viewer webserver [73], KEGG (Kyoto Encyclopedia of Genes and Genomes) Automatic Annotation System (KAAS) database [74] and the bacterial version of the antiSMASH webserver [75]. All options were selected with the default parameters. Additionally, for deeper metabolic investigation, the amino acid sequences were subjected to analysis with BLAST-KOALA [76]. KO (KEGG Orthology) assignments were performed using a modified version of the internally used KOALA (KEGG Orthology And Links Annotation) algorithm (BLASTKOALA), after a BLAST search against a non-redundant dataset of pangenome sequences [76]. Hits showing at least 50% identity with the reference protein were considered significant. Each hit was verified manually using BLASTp analysis.
The analysis of the spread and conservation of proteins encoded within the core of identified fosmids was performed throughout the protein search of proteins from within the core against the protein non-redundant (nr) National Center for Biotechnology Information (NCBI) database (accessed on April 14, 2022) with the application of diamond v2.0.13 using the following parameters: -ultra-sensitive -query-cover 75 -subject-cover 75 -id 75 -max-target-seqs 0 [77]. Based on the subject proteins' accession numbers, the location of their genes was accessed with NCBI E-utilities. Only INSDC databaseoriginating records were considered. The alignment of selected contigs was performed with clinker v0.0.20 with a 75% sequence identity threshold [78].

Fosmid library preparation and DNA cloning
The ANT_H4 DNA of 30-50 kb was end-repaired to blunt end and cloned into the copy control pMPO579 vector using the CopyControl Fosmid Library Production Kit, according to the manufacturer's protocol [41,79].
For subcloning, amplification of the relevant DNA region of the previously isolated ANT_H4-derived fosmid was performed by PCR using synthetic oligonucleotide primers, dNTPs and HiFi polymerase (Qiagen; with supplied buffer) in a Mastercycler (Eppendorf ). The forward and reverse primers used in this study were 5′CGT ACT GCA GGT CGT TAC GGT TCA TCT TGT G and 5′GAC TCT CGA GTG TAA TAG GCC GTT ACC AGT C, respectively. The PCR product was then analyzed by electrophoresis on 0.8% agarose gel. The confirmed product of 1,689 bp was then digested with XhoI and PstI restriction enzymes and cloned into the multiple cloning site located within the lacZα sequence of pBluescript II SK to facilitate selection by standard blue-white screening. Further transformation into E. coli TG1 was performed according to the Kushner method [62].

Quantification of melanin-like compounds
The ANT_H4 strain was cultivated overnight in lysogeny broth (LB) medium at 18 °C with rotary shaking set to 3 RCF. Cultures were then centrifuged (4000 RCF for 5 min) and washed three times with 0.85% saline solution. Next, bacteria were added in triplicate to the fresh LB medium supplemented with 0.1% (w/v) L-tyrosine and cultivated at 18 °C with rotary shaking set to 3 RCF. The initial optical density at 600 nm (OD 600 ) was 0.1. The idiophase was established after 5 days of cultivation (confirmed by the number of colony-forming units (CFU) counted after two consecutive days) and the experiment was continued up to day 10 to maximize metabolite production. Bacteria were then centrifuged (4000 RCF for 5 min) and the supernatant was measured spectrophotometrically at a wavelength of 400 nm [30]. A standard curve was prepared with polymerized homogentisic acid (Sigma-Aldrich) while the LB medium was a blank. Purification of the metabolite followed from the bacterial-free supernatant which was acidified with HCl solution to pH 2 and left overnight for melanin precipitation. In the next step, the residue was washed three times with ddH 2 O and ethanol. To reduce runaway esterification reactions, the metabolite was dried out of the alcohol. To dissolve the compound, the precipitate was suspended in ddH 2 O at pH 11 (NaOH solution). For the majority of experiments, the pH was adjusted to 7, which did not affect the solubility of the compounds. In order to remove the salt, the resuspended metabolite was then dialyzed in 10 kDa dialysis tubing for 24 h against water [15,30].

High-performance liquid chromatography (HPLC)
The purified and water-dissolved melanin-like compound from ANT_H4 was screened using an Agilent 1100 HPLC system (Agilent Technologies) with a DAD UV detector. Multiple wavelength monitoring was performed at λ = 216, 254, and 310 nm. 10 µl of the sample was injected onto an HPLC column (100 × 4.60 mm) filled with 5 µm Luna-100 C-18 (Phenomenex, Aschaffenburg, Germany). The samples were analyzed by linear gradient elution using H 2 O + 0.1% formic acid as solvent A and acetonitrile + 0.1% formic acid as solvent B at a flow rate of 2 mL/min. For solvent B, the gradient was from 5 to 100% in 10 min with a 3 min isocratic elution at 100% [80].

Fourier-transform infrared (FT-IR) spectroscopy
Infrared spectra were recorded using the Nicolet iS50 FT-IR spectrometer (Thermo Scientific, Waltham, MA, United States) with a DTGS detector. Sample preparation was carried out as follows: 2 mg of dried melanin compound was mixed with 230 mg of KBr (Sigma-Aldrich) and carefully ground with agate mortar. Pellets were prepared using a press under a pressure of 10 T for 20 s. To collect spectra, the transmission mode was used with a resolution of 4 cm −1 and typically 256 scans were taken for each sample. A KBr pellet was used as a reference background [81].

Free radical scavenging activity
The DPPH (2,2-diphenyl-1-picrylhydrazyl) method was used to measure the free radical scavenging activity. The dilutions were prepared as follows: 2 mL of 0.1 mM DPPH in methanol was added to 2 mL of methanol containing different amounts of melanin-like compound, i.e., to reach final concentrations of 1, 0.5, 0.25, 0.1, 0.01, and 0.001 mg/mL. The absorbance at 517 nm was measured spectrophotometrically (Evolution 260 Bio, Thermo Fisher Scientific) after 30 min. The scavenging of the DPPH radical (%) was calculated according to the formula ((A0 − A1)/A0 × 100) [82], where A0 is the absorbance of the control reaction and A1 is the absorbance of reactions containing melanin from the ANT_H4 strain. Each experiment was performed in triplicate.

Sun protection factor (SPF) analysis
0.1 g of commercial sunscreen marked with an SPF of 5, 15, and 25, respectively, was added to 10 mL of absolute ethanol. To this solution, a melanin-like compound was added to final concentrations of 0.1 and 0.01 mg/mL. A variant without supplementation was used as a control. The absorbance of the mixture in the UV range (290-320 nm) was quantified at 5 nm intervals using ethanol as the blank. SPFs were calculated according to the Mansur method [83,84]. Each experiment was performed in triplicate.

Bioconsolidation assay
In order to establish the ability of the melanin-like compound to increase iron retention, 10 mL of two water solutions of melanin (0.1 and 0.01 mg/mL) were incubated on a mechanical cradle with 1 g of hematite (Geogut, Poland) and pyrite (Geogut) for 1, 24, and 48 h. The extracts were then filtered through 0.22-µm filters and submitted for further analysis. Distilled water was used as a control extraction solution. All extractions were performed in triplicate. The initial pH of the melanin solutions was 7.

Analysis of iron concentration in the samples
The amount of iron was measured by Graphite Furnace Atomic Absorption Spectroscopy (GFAAS) using a Thermo Scientific SOLAAR M Series (TJA Solution, SOLAAR M, Cambridge, United Kingdom). The gas mixture was air and acetylene. The graphite cuvette duty cycle was as follows: evaporation-100 °C/30 s; incineration-1100 °C/20 s; atomization2,2-diphenyl-1-picrylhydrazyl 2100 °C/3 s; cleaning-2500 °C/3 s, and the calibration curve range was 0-20 µg/l with a lower limit of quantification of 0.1 µg/l. A deuterium lamp (TJA Solution) was used for background correction. Iron standard solutions (Merck, Darmstadt, Germany) were prepared in 3% HNO 3 .

Iron chelation assay
In order to establish chelating properties, 0.01, 0.1, 0.25, 0.5, 1, 2, and 4 mg/mL of water solutions of pyomelanin were added in a 1:1 ratio to the CAS (chrome azurol S) reagent [66]. Distilled water was used as a negative control. The standard curve (R 2 = 0.99) was prepared against deferoxamine mesylate salt (Sigma-Aldrich). All experiments were performed in triplicate. After an hour of incubation, the absorbance at 630 nm was measured using an automated microplate reader (Sunrise TECAN, Tecan Trading AG, Männedorf, Switzerland).

In vitro hairy roots cultivation
Calendula officinalis L. hairy root cultures were produced in a ½ MS medium at a temperature of 24 °C, in the dark, on a shaker (1.5 RCF). Weighed samples of 3-week-old roots were transferred to a fresh medium, to which an aqueous solution of pyomelanin was added at a concentration of 0.01% (w/v). The control variant was distilled water. The tests were carried out in triplicate per variant after: (T1) transferring to a new medium with an elicitor, (T2) 7, and (T3) 14 days.

Extraction of primary and secondary metabolites from hairy roots cultures
The hairy roots were weighed and extracted with: (i) diethyl ether for the determination of free sterols and oleanolic acid (pentacyclic triterpenoid) and (ii) methanol to investigate sterol glycosides and oleanolic acid saponins. The post-culture medium was also tested for the presence of isolated oleanolic acid saponins. The extraction of plant metabolites from hairy roots tissues was carried out in Soxhlet apparatuses. Saponins (extracted with butanol from post-culture medium as well as saponins from hairy root tissue) and sterol glycosides underwent acidic hydrolysis to release aglycones [60]. They were then fractionated against sitosterol and oleanolic acid (Sigma-Aldrich) standards using thin layer chromatography (TLC) in an appropriate solvent system (chloroform: methanol 97:3 v/v). Prior to analysis by gas chromatography with mass spectrometry (GC-MS), fractions containing oleanolic acid were additionally methylated [59].

Identification and quantification of triterpenoids by gas chromatography coupled to a flame ionization detector and mass spectrometer (GC-FID/MS)
An Agilent Technologies 7890 A gas chromatograph equipped with a 5975C mass spectrometric detector was used for qualitative and quantitative analyses. Samples dissolved in diethyl ether:methanol (5:1, v/v) were applied (in a volume of 1-4 μL) using a 1:10 split injection. The column used was a 30 m × 0.25 mm (L × I.D), 0.25 μm particle size (HP-5MS UI, Agilent Technologies, Santa Clara, CA, USA). Helium was used as the carrier gas at a flow rate of 1 mL/min. The separation was made either under isothermal conditions at 280 °C or at the programmed temperature; an initial temperature of 160 °C held for 2 min, then increased to 280 °C at 5 °C/1 min, and the final temperature of 280 °C held for a further 44 min. Additional parameters were employed as follows: inlet and FID temperature 290 °C; MS transfer line temperature 275 °C; quadrupole temperature 150 °C; ion source temperature 230 °C; EI 70 eV; m/z range 33-500; FID gas (H2) flow 30 mL/min (hydrogen generator); and airflow 400 mL/min. Individual compounds were identified by comparing their mass spectra with library data from Wiley 9th ED. and NIST 2008 Lib. SW Version 2010 or previously reported data, and by comparison of their retention times and corresponding mass spectra with those of authentic standards, when available. Quantitation was performed using an external standard method based on calibration curves determined for the compounds belonging to representative triterpenoid classes; α-amyrin for triterpene alcohols, oleanolic acid methyl ester for triterpene acid methyl esters, and sitosterol for steroids.

Statistical analysis
The significance of the differences between the mean values of the control and treated plant and soil samples was statistically evaluated by a two-tailed t-test at p ≤ 0.05. The Mann-Whitney-Wilcoxon test was applied whenever data failed to present a normal distribution or had different variances. The statistical analysis was carried out using the XLSTAT (version 2022.1) program.

Supplementary Information
The online version contains supplementary material available at https:// doi. org/ 10. 1186/ s12934-022-01990-3. Table S1. GC-MS data (retention times and characteristic ions of mass spectra) of identified steroids and triterpenoids. Table S2. The content of steroids and triterpenoid alcohols in hairy roots tissue. Table S3. The content of sterol glycosides in hairy roots tissue. Table S4. Free oleanolic acid (OA) content in hairy roots tissue. Table S5. Oleanolic acid saponins (OA) released to the medium. Table S6. Oleanolic acid saponins (OA) content in hairy roots tissue.